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Abstract. We reanalyze the inference of the protosolar 
abundance of deuterium made by Geiss (1993) from mea- 
surements of ( 3 He/ 4 He)sw m the solar wind. We use an 
evolutionary solar model with microscopic diffusion, con- 
strained to fit the present age, radius and luminosity, as 
well as the observed ratio of heavy elements to hydrogen. 
The protosolar ( 2 H/ 1 H) P is obtained from the best fit of 
( 3 He/ 4 He)sw- Taking for the protosolar ( 3 He/ 4 He) p the 
value measured in Jupiter by the Galileo probe (Niemann 
et al. |l996l), we derive ( 2 H/ 1 H) P = (3.01 ± 0.17) x 10~ 5 . 



Compared to the present interstellar medium value (Lin- 
sky et al. 1993), this result is compatible with models of 
the chemical evolution of the Galaxy in the solar neigh- 
borhood; it is also marginally compatible with the Jovian 
( 2 H/ 1 H) J = (5 ± 2) x 10~ 5 ratio measured by Galileo. 



Key words: Nuclear reactions, nucleosynthesis, abun- 
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1. Introduction 



Since the pioneering work of Geiss & Reeves ( 1972| ), a 
number of papers have been devoted to the fundamental 
question of the estimate of the proto solar abundance of 
deuterium as reviewed by Geiss ( 1993 ), hereafter referred 
as G93. The current value of ( 2 H/ 1 H) P = (2.6±1.0) x 10~ 5 , 
is that derived by G93 from a reanalysis of the measure- 
ments of ( 3 He/ 4 He) SW in the Solar Wind (SW). 

At the end of the pre-main sequence the young Sun is 
still fully convective with a temperature, at center, of a 
few million K; the mixing makes that all the deuterium 
is converted to 3 He via the reaction 2 H(p, 7) 3 He. As 



pointed out by Geiss k Reeves (1972), if ( 4 He/ 1 H)= 0.1 



is no longer changed at the surface then, ( 2 H/ 1 H) p is 
equal (G93) to ( 4 He/ 1 H) p times the difference between 



( 3 He/ 4 He) Q presently observed at the surface and the pro- 
tosolar ( 3 He/ 4 He) p : 



2H 



3 He\ / 3 Hc\ 



( 3 He/ 4 He) q is measured only in the solar wind, from the 
data available G93 concluded to ( 3 He/ 4 He)sw= (4-5 ± 
0.4) x 10~ 4 . However G93 analyzing the various processes 
which may fractionate the helium isotopes in the interior 
and in the atmosphere of the Sun, has estimated that the 
actual ratio to be considered for inferring ( 2 H/ 1 H) p is the 
solar wind value divided by a factor 1.1 ± 0.2, resulting 
in ( 3 He/ 4 He) = (4.1 ± 1.0) x lO" 4 . ( 3 He/ 4 He) p has been 
estimated by G93 from meteoritic data to be equal to 
( 3 He/ 4 He) p = 1.5 x 10~ 4 . 

An other way to determine the protosolar deuterium 
abundance is to evaluate ( 2 H/ 1 H) i y in Jupiter. This ap- 
proach is based on the fact that Jupiter is mainly made of 
hydrogen which originates from the primordial solar neb- 
ula. According to current models of formation of Jupiter, 
some amount of ices, presumably enriched in deuterium 
with respect to the protosolar value, may have been mixed 
to hydrogen during the planetary formation, but their 
mass is currently considered as too small to have sig- 
nificantly incr eased ( 2 H/ 1 H) t7 (Hubbard & MacFarlane 



1980). Although the various estimates of ( 2 H/ 4 H)j- from 
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remote sensing observations of HD and of CH3H arc un- 
certain, it was generally considered, prior to the arrival of 
the Galileo mission to Jupiter, that the Jovian deuterium 
abundance is consistent with the G93 value (Lecluse et 
al. 1996| ). The preliminary result of the mass spectrome- 
ter aboard the Galileo atmos pheri c probe - ( 2 H/ 1 H)j-= 
(5 ±2) x 10~ 5 (Niemann et al. |l996j) - is however substan- 
tially higher than the G93 value. This result, if confirmed, 
would have profound cosmological implications. 

Compared to the present day InterStellar Medium 
(ISM) value of ( 2 H/ 1 H)i MS = (1-47 - 1.72) x lO" 5 (Lin- 
sky et al. 1993 ), it suggests a deuterium destruction in 



2 



D. Gautier & P. Morel: A Reestimate of the Protosolar ( 2 H/ 1 H) ratio 



4.55 Gyr, the solar age, inconsistent with current models 
of chemical evolution of galaxies. Alternatively, it might 
be that ( 2 H/ 1 H)ims varies upon the line of sight so that, 
the concept of a well defined ISM deuterium abundance is 
not reliable, as advocated by some authors e.g., Ferley et 
al. ( 1995D - A third possibility is that the current protosolar 
value ( z H/ 1 K)p, is underestimated. 

Two circumstances prompted us to reanalyze the de- 
termination of G93: - First, Niemann et al. ( 1996| ) have 
measured in Jupiter ( 3 He/ 4 He)j= 1.1 x 10~ 4 which must 
be the exact ratio in the primitive solar nebula. As pointed 
out by these authors, this results in an increase of the es- 
timation of ( 2 H/ 1 H) p . - Second, since 1993 the introduc- 
tion of microscopic diffusion in solar modeling have signif- 
icantly increased the accuracy of models (e.g., Basu et al. 
199§). A consequence of the microscopic diffusion is the 



strumcnt aboard the Ulysses spacecraft (Bodmer et al. 



change, with respect to time, of surface values of ( 4 He/ 1 H) 
and ( 3 Hc/ 4 He). A result directly relevant to this Letter is 
the inferred helium mass fraction at the surface namely 
Y Q = 0.2 42 ± .003 (Perez Hernandez & Christensen- 
Dalsgaard 1994 ) substantially less than t he p rotosolar 
value Y p = 0.275 (Bahcall & Pinsonncault |1995| ). There- 
fore, we propose here to complete the approach of G93 
by using a solar evolutionary model with microscopic dif- 
fusion. The modeling of the solar evolution and the data 
employed are described in Section ||; Section || is devoted 
to results and discussions. We conclude in Section |]. 

2. Solar modeling and data 

The calibrated solar models discussed in this paper in- 
clude the pre-main sequence evolution. The calibration 
consists in the adjustment of, i) the ratio of the mixing- 
length to the pressure scale height, ii) the initial mass frac- 
tion Xp of hydrogen and, iii) the initial protosolar mass 
fraction of heavy elements to hydrogen, in order that the 
models have at present solar age, the observed luminosity, 
radius (Guenter et al. 1992Q and mass fraction of heavy 
elements to hydrogen (Grevesse & Noels 1993). The mod- 



els h ave been computed using the code CESAM (Morel 
1997 ). The relevant features are: the changes due to nu- 
clear reactions, microscopic diffusion and convective mix- 
ing are explicitly computed for 1 H, 2 H, 3 He, 4 He, 7 Li, 
7 Be, 12 C, 13 C, 14 N, 15 N, 16 and 17 which enter into 
the most important nuclear reactions of the PP+CNO 
cycles; the protosolar abundance of each heavy element 
is derived from X p according to the nuclide abundances 
of Anders & Grevesse (1989); ( 3 He/ 4 He) p is taken equal 



to ( He/ He) j— 1.1 x 10 as measured by the Galileo 
probe into Jupiter; the microscopic diffusion coefficients 
of Michaud & Proffitt fll993| ) are used, all species, but 
1 H and 4 He, are trace elements. ( 2 H/ 1 H) p is taken as a 
parameter and constrained to fit ( 3 He/ 4 He) . 

As pointed out by G93, 3 He could be somewhat favored 
over 4 Hc in the solar wind acceleration process (Burgi & 
Geiss 1986 ). However, the first results of the SWICS in- 



1995) do not provide evidence that such a fractionation 



could exceed the uncertainties on the measurements of 
the helium isotopic ratio. Moreover, the chromospheric 
process is not expected to result in a significant mass frac- 
tionation (G93). Therefore we assume that ( 3 He/ 4 He) is 
equal, within uncertainties, to ( 3 He/ 4 He)sw- The values 



Fig. 1. Profiles of ( 3 He/ 4 He) in envelopes of solar models com- 
puted with ( 2 H/ 1 H) P : 3.2 10~ 5 (thick), 3.6 10" 5 (medium) and 
2.8 10 -5 (thin). The observed ( 3 He/ 4 He) S w from APOLLO 
(1), ISEEE-3 (2) and Ulysses (3) are plotted with their error 
bars. 

of ( 3 He/ 4 He)sw obtained by three experiments on three 
different space missions are: 

(1) Apollo ( 3 He/ 4 He) sw = (4.25 ±0.21) x 10~ 4 (Geiss et 



al. 1972) 



(2) ISSEE-3 ( 3 Hc/ 4 Hc) sw = (4.88 ±0.48) x 10" 4 (Coplan 



et al. 1984, Bochsler 1984), 



(3) Ulysses ( 3 He/ 4 He) 



sw= 



et al. 1995) 



(4.4 ± 0.4) x 10" 4 (Bodmer 



3. Results and Discussion 

For models initialized for various values of ( 2 H/ 4 H) p , the 
( 3 He/ 4 He) profiles are plotted Fig. | for R £ O.55i? and 
compared to the three observed ( 3 He/ 4 He)sw values. Ac- 
cording to these models, ( 2 H/ 4 H) p can be fitted by the 
quadratic polynomial: 

( 2 H/ 1 H) p = 2.675 x 10~ 5 + (9.364 x 10~ 2 - 2.452a;)a;. (1) 

here x = ( 3 He/ 4 He) - 4.0 x 10~ 4 , ( 3 He/ 4 He) profiles 
within < R < 1.2R Q are plotted Fig. | at various ages 
of a solar model computed with ( 2 H/ 1 H) p =3.01 x 10~ 5 . 
During the pre-main sequence 2 H is converted to 3 He and, 
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Fig. 2. ( 3 He/ 4 He) profiles for R < 1.2R Q , in a solar 
model initialized with ( 2 H/ 1 H) P = 3.01 x 1CT 5 for various 
ages: OMyr (dash-dot-dot-dot), 0.138 Myr (dashed), 20Myr 
(dash-dot-dash), 49.1 Myr (thin full), 200 Myr (dotted) and 
4550 Myr (thick full). 



Fig. 3. Changes of ( 3 He/ 4 He) (thick) and 4 He/ 1 H (thin) ra- 
tios (with different scaling) at surface with respect to time for 
the solar model of Fig. ^. The dashed heavy line (49 Myr) sep- 
arates the pre-main sequence (PMS) from the main sequence 
(MS). 



due to mixing, ( 3 He/ 4 He) increases through the model. 
At time t =~ 20 Myr the convection zone has receded 
almost to its present days location, at center, the tem- 
perature is not high enough to convert 3 He into 4 He via 
He 3 (He 3 ,2p) 4 He, then ( 3 He/ 4 He) is maximum there. At 
t =~ 49 Myr i.e., zero age main sequence, the flat profile 
of ( 3 He/ 4 He) for radius R ig O.li?0, is due to the mix- 
ing in the convective core resulting form the conversion 
of 12 C into 14 N; then the nuclear reactions reach their 
equilibrium. Around the center the increase of 4 He de- 
presses ( 3 He/ 4 He) and its maximum progressively reaches 
its present days location around 0.3i? Q ; the gravitational 
settling being more efficient for He than for 3 He, at sur- 
face ( 3 He/ 4 He) slowly increases until present days. As 
seen in Fig. along the evolution, ( 3 He/ 4 He) varies from 
1.10 x 10~ 4 to 4.361 x 10~ 4 and ^He/^H) from 0.0966 to 
0.0830. 

Using Eq. [l] the three observed ( 3 He/ 4 He)sw values 
allow to infer: 

(1) ( 2 H/ 1 H) p = (2.91 ±0.19) x 10- 5 , 

(2) ( 2 H/ 1 H) p = (3.50 ±0.43) x 10" 5 , 

(3) ( 2 H/ 1 H) P = (3.05 ±0.36) x 10~ 5 . 

Assuming no systematic errors, the weighted mean of 
these determinations results in: 

( 2 H/ 1 H) p = (3.01 ± 0.17) x 10" 5 . 



3 He/ 4 He) p = 0.0966 value, they would have found 



Niemann et al. (|1996| ) have derived a similar value but us- 
ing our 

( 2 H/ 1 H) p = (2.9 ±1) x 10 

Our new estimate of ( 2 H/ 1 H) p overlaps within the do- 
main of uncertainty the Galileo result. At this point, we 



can thus consider that it is consistent with Galileo. How- 
ever, it rules out Jovian values higher than ( 2 H/ 1 H)j-^ 
3.18 10 -5 if we assume that Jupiter is representative of 
the isotopic composition of the nebula. As mentioned in 
Sect. |H ( 2 H/ 1 H) J results from a mixing of hydrogen origi- 
nating from the nebula - and thus in protosolar abundance 
- with ices more or less enriched in deuterium. Assuming 
that the two reservoirs equilibrated at high temperature 
at the time of the formation of the planet, the present 
hydrogen in Jupiter may have been somewhat enhanced 
in deuterium if the amount of ices was large enough. This 
question deserves to be reexamined in the light of the most 
recent models of interiors of Jupiter. Following Hubbard 
& MacFarlane (1980), the deuterium enhancement is a 
function of the mass Mi of the ices embedded in Jupiter 
and of their deuterium enrichment / with respect to the 
protosolar value. According to Guillot et al. (1997), the 
amount of ices should not exceed 32 Earth masses in the 
extreme case. The value of the enrichment / depends on 
the origin of Jovian protoices. They may have been formed 
in the solar nebula, which implies that / not to exceed 2.5. 
In such a case, the deuterium enhancement in Jupiter is 
negligible. Alternatively, ices may have originated directly 
from the protosol ar clou d following the scenario discussed 
by Lunine et al. ( 1991 ) and would have then kept their 
interstellar isotopic signature. The recent determin ations 
of ( 2 H/ 1 H) in water in P/Halley (Eberhardt et al. |l995|) 
and in Hyakutake comet (Gautier et al. 1996| ) suggest / 
of the order of 10. This maximum scenario then results in 
a deuterium enhancement of 20%. Accordingly, ( 2 H/ 1 H)j 



4 



D. Gautier & P. Morel: A Reestimate of the Protosolar ( 2 H/ 1 H) ratio 



should then be equal to ( 2 R/ 1 H) J = (3.65 ± 0.2) x 1(T 5 , 
entirely within the error bars of the Galileo value. 

Ground based remote sensing determinations of 
( 2 H/ 1 H)j- conclude to values less than 3.0 x 10~ 5 (Lecluse 
et al. 1996D , in conflict with the maximum enrichment 
scenario. The preliminary analysis of observat ions from 
ISO suggest similar results (Encrenaz et al. 1996). In 



contrast, the reanalysis o f the Voyager infrared observa- 
tions by Carlson et al. ( 1993j ) results in ( 2 H/ 1 H)j_be- 
tween 3.7 x 10~ 5 and 6.0 x 10~ 5 (Lecluse et al. |l996| ), 
a value compatible with the maximum enrichment sce- 
nario. Preliminary results of HST observations of th e Ly- 
man a emission of Jupiter (BenJaffel et al. 1996 ) give 
( 2 H/ 1 H) J =(5.9 ± 1.4) x 10~ 5 , a value much higher than 
the upper limit of the maximum enrichment case. 

In summary, we cannot yet decide whether the abun- 
dance of deuterium in Jupiter is representative of the pro- 
tosolar value or if it has been somewhat enhanced during 
the planetary formation. 

The lower limit of our revised protosolar ( 2 H/ 1 H) p va- 
lue, namely ( 2 H/ 1 H) p = 2.83 x 10~ 5 , is higher than the up- 
per limit of the ISM determination of Linsky et al. ( |1993 ), 
by a factor 1.6. This is consistent with the evolution in 
4.55 Gyr of the deuterium abundance in the solar neigh- 
borhood, estimated from models of chemical evolution of 
the Galaxy with infall of primordial composition (Prant- 
zos 1996). Obtaining higher depletion factors requires dif- 



ferent models which invo ke ga lactic winds in the Galaxy 
(Vangioni-Flam & Casse 1995 ). 



4. Conclusion 

We propose for the protosolar deuterium abundance: 
( 2 H/ 1 H) p = (3.01 ± 0.17) x 10~ 5 . This result is based 
on observed values of ( 3 He/ 4 He)sw m the solar wind and 
of ( 2 H/ 1 H)j- on Jupiter; it takes into account the micro- 
scopic diffusion of elements in solar modeling. Systematic 
errors could occurs: i) although not yet detected, a frac- 
tionation between 3 He and 4 He could occur between the 
surface and the solar wind, resulting in a small decrease of 
the inferred ( 2 H/ 1 H) p , ii) microscopic diffusion coefficients 
could be revised, affecting the calculation of ( 3 He/ 4 He)0. 
The new protosolar value is consistent with the lower 
range of the preliminary determination of ( 2 H/ 1 H)j- in 
Jupiter obtained by the Galileo mass spectrometer. The 
difference between our value and the present ISM value 
of Linsky et al. ( |1993[ ) is compatible with some models 
of the chemical evolution of the Galaxy. Our result does 
not significantly changes the sum ( 2 H/ 1 H) p +( 3 He/ 1 H) p = 
4.07 x 10~ 5 , from its previous estimate of 4.1 x 10 -5 by 
G93. In fact, using our value ( 4 He/ 1 H) p = 0.0966, this au- 
thor would have found 4.05 x 10~ 5 . 
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